Mass Limits on Neutralino Dark Matter 
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We set an upper limit on the mass of a supersymmetric neutralino dark matter particle using the 
MicrOMEGAS and DarkSUSY software packages and the most recent constraints on relic density 
from combined WMAP and SDSS data. We explore several different possible scenarios within the 
MSSM, including coannihilation with charginos and sfermions and annihilation through a massive 
Higgs resonance, using low energy mass inputs. We find that no coannihilation scenario is consistent 
with dark matter in observed abundance with a mass greater than 2.5 TeV for a wino-type particle 
or 1.8 TeV for a Higgsino-type. Contrived scenarios involving Ifiggs resonances with finely-tuned 
mass parameters can allow masses as high as 34 TeV. The resulting gamma-ray energy distribution 
is not in agreement with the recent multi-TeV gamma ray spectrum observed by H.E.S.S. originating 
from the center of the Milky Way. Our results are relevent only for dark matter densities resulting 
from a thermal origin. 



I. INTRODUCTION 

The existence of cold dark matter in the universe 
is well-established by recent astrophysical observations. 
While the particle nature of dark matter remains a mys- 
tery, its effect of promoting structure formation via hier- 
archical formation of gravitational potential wells from 
primordial density fluctuations is well-documented in 
both high resolution N-body simulations and observa- 
tions of large-scale structure [11, . The properties re- 
quired for a thermally-produced particulate cold dark 
matter include high mass, lack of strong or electromag- 
netic couplings, and stability over the lifetime of the uni- 
verse. No standard model particle is capable of simul- 
taneously satisfying these demands. Minimal supersym- 
metry provides a number of viable candidates among its 
spectrum of fundamental particles 3, 4]. The lightest of 
the four neutralinos, which in models we will consider is 
also the lightest supersymmetric partner or LSP, consists 
of various linear combinations of the neutral bino, wino, 
and Higgsino states. The LSP has been proposed as an 
excellent candidate not only because of the above prop- 
erties, but also because the expected cross section gives 
rise naturally to the observed dark matter mass density 
during thermal decoupling. If supersymmetry is to solve 
the gauge hierarchy problem, it should be broken in such 
a way that the partner particles attain mass corrections 
on the order of the electroweak scale, ~ 100 GeV. A 
weakly-interacting dark matter particle has the correct 
cross section to produce the observed relic density. 

Attempts to discover the particle nature of dark mat- 
ter take on a three-pronged approach. Terrestrial ex- 
periments to directly detect dark matter particles pass- 
ing through the earth are underway [1, [a, 0| > ^'^'^ it is 
expected that the next generation of particle accelera- 
tors will be capable of producing a weakly interacting 
particle with a mass near the electroweak scale Q. In 
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contrast to direct detection and production of dark mat- 
ter, indirect detection techniques search for the products 
of dark matter annihilation. Self-annihilations between 
clustered particles are expected to produce a variety of 
high-energy cosmic rays, photons, and neutrinos. Sites 
for these events include the central den sity spikes of dark 
matter halos [1, [l^, [H [H, [ll H [H, [11 , diffuse radia- 
tion from the halo at large [iTI . |l8| , substructure within 
the galactic halo [l^ [l^l , and within astronomical bodies 
such as stars and planets, includin g ou r own earth and 
sun [U, [23, [2^ . As mentioned in |l3j , the gamma ray 
flux per energy bin from neutralino annihilations in a re- 
gion of space requires inputs from two factors, one being 
the cross section times the expectation value for the num- 
ber of gammas produced, and the other an astrophysical 
variable determined by the distance and density proflle 
of the region. In addition to assuming a supersymmetric 
model, calculation of the signal from the center of our 
galaxy requires knowledge about the proflle and normal- 
ization of the inner halo. On scales larger than about 1 
kpc, it is well-established through N-body simulations 
that the halo has a p ower-law radial density profile with 
index -1 to -1.5 [24 [25l [26l Frequent gravitational 
scattering events between dark matter and stars leads to 
an equilibrium profile of p ^ 7-^3/2 jj^ the inner 2 parsecs 
where dark matter is a negligible proportion of the 
total mass. Since the normalization is indeterminate, as 
a fiducial model we interpolate inwards from an NFW 
proflle [IJ]. While this may be a reasonable assumption, 
it is likely that baryonic compression, in which baryonic 
matter losing energy through radiative processes falls in- 
ward and consequently redistributes dark matter, has in- 
creased this number, possibly by several orders of mag- 
nitude [Tsj . For a weakly-interacting particle of mass 
greater than 100 GeV, it is plausible that a dark matter 
signal will be observed from the Galactic center by the 
current generation of ground-based gamma ray detection 
experiments, provided that baryonic compression has in- 
creased the central density by at least a modest amount 
over the flducial value [l^l, although astrophysical pro- 
cesses could be a potentially troubling background for 
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these searches 

A recent observation of a TeV gamma-ray signal 
from the Galactic center by the H.E.S.S. atmospheric 
Cherenkov telescope [1^, [lo] has motivated this deter- 
mination of a theoretical upper limit on the mass of 
neutralino dark matter. This result was confirmed re- 
cently by the MAGIC telescope [31]. The signal was also 
observed previously by the CANGAROO-II experiment 
[s^ . though this result was inconsistent with the newer 
H.E.S.S. results; the large flux observed by CANGAROO 
at low energies was not seen by H.E.S.S. All groups have 
stated that the signal is consistent with a point source, 
and H.E.S.S. and MAGIC agree on a logarithmic slope of 
about -2.25. Without a clear indication of an annihilat- 
ing particle, such as line features in the spectrum, or an 
observation that disfavors an annihilation scenario, such 
as time variability, the source of the signal remains un- 
certain, although the extended power law nature of the 
observed spectrum does not fit well with the expected 
rollover shape of an annihilation spectrum. An analysis 
by the H.E.S.S. group of their 2004 data, which extended 
to approximately 30 TeV, was unable to find any an- 
nihilation spectra which reproduced the observed power 
law, and they proposed that the signal must be primarily 
non-dark matter in origin. While the newest data from 
H.E.S.S. does not seem to be consistent with the annihi- 
lation spectrum of either a supersymmetric neutralino or 
Kaluza-Klein dark matter particle [33| , models in which 
a signal from annihilating dark matter is masked by emis- 
sions of astrophysical sources are still possible [30j. An- 
other analysis of the 2004 data by Profumo searched for 
optimal spectral fits based on final state channels; it was 
determined that the dark matter annihilation remains 
a possible interpretation of the H.E.S.S. data for a re- 
stricted set of final states Qg. The 2003 H.E.S.S. data, 
which extends to 9 TeV, can similarly be fitted with fewer 
constraints on the final states of the annihilating parti- 
cle. Mambrini and collaborators [3^ searched for a neu- 
tralino annihilation spectrum in a non-universal super- 
gravity model which could fit the 2003 data. They were 
successful in finding reasonable fits, though none of the 
points in parameter space for their high mass candidates 
were consistent with WMAP constraints. Possible astro- 
physical sources of TeV gamma-rays include jets or the 
shocks in the accretion flow into the central black hole 
[36l [srj . Another possibility, which may be ruled out by 
lack of time variability, is a signal from particles acceler- 
ating near the event horizon of a rotating super-massive 
central black hole [2i,|3i|. 

While various authors have conducted surveys of su- 
persymmetric parameter space while categorizing dark 
matter candidates (e.g. [H, HO, Ej), few papers have 
explicitly searched for upper mass limits on these candi- 
dates. The effects of coannihilation in a low-energy ef- 
fective MSSM model were considered in [i^, which cal- 
culated the degree to which relic density was reduced 
by various channels. We conduct a similar survey, but 
with the goal of examining the regions of maximal mass 



in Higgsino coannihilation. In [i^, the DarkSUSY soft- 
ware package was used to explore coannihilation in the 
mSUGRA parameter space. In the chargino coannihila- 
tion region, the LSP is mostly Higgsino, and coannihi- 
lates with a nearly-degenerate chargino and second Hig- 
gsino. This publication reported a cosmological limit 
of ^1500 GeV resulting from chargino coannihilation, 
which is nearly consistent with our findings for a pure 
Higgsino in this region. They also examined coannihi- 
lations of a bino-type LSP with sleptons and put forth 
the claim that coannihilation processes do not allow ar- 
bitrarily high masses, in contrast to some previous au- 
thors. Higgsino dark matter in an mSUGRA framework 
was also considered in , where a WM AP-favored mass 
range of approximately 1 TeV is found. Wino-type dark 
matter appears in minimal anomaly-mediated supersym- 
metry breaking (mAMSB) models, and has a cosmologi- 
cal mass bound of over 2 TeV . Profumo [33| looked 
at several different scenarios, including a mAMSB wino, 
annihilation through a heavy Higgs resonance, and QCD 
effects in gluino annihilations which can allow dark mat- 
ter masses in the hundreds of TeV. 

We use a model-independent approach in the calcu- 
lation of the rriLSP upper bound, one which could be 
applicable to many individual supersymmetry-breaking 
models. To this end, we have used model inputs at the 
electroweak scale, allowing us to control individual in- 
puts to the mass parameters and compute the dark mat- 
ter relic density ^Ilsp- Calculations were done using the 
MicrOMEGAs [46] and DarkSUSY [47,] software pack- 
ages. Both of these codes compute SUSY relic density 
via a numerical solution of the Boltzmann equation, in- 
cluding the cross sections of any relevant coannihilation 
processes. Both programs also provide several options 
for inputting the supersymmetric particle spectrum, in- 
cluding electroweak and GUT-scale minimal inputs, or 
determining the mass of each particle independently. Mi- 
crOMEGAs calculates the contribution of individual de- 
cay channels to standard model products, as a function of 
the contribution of each channel to , there being an 
approximate inverse relationship between the total cross 
section and final relic density. DarkSUSY also provides 
the necessary tools to calculate the flux of high energy 
gamma-rays from halo annihilation, as well as a variety 
of other direct and indirect detection signals. 



A. The Boltzmann Equation 

The density evolution of any particle x iu the thermal 
bath of the early universe is governed by the Boltzmann 
equation: 




Here a is the cosmological scale factor, and n and n^"^ 
are the number density and equilibrium number density 
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of the particle species. The thermally-averaged cross sec- 
tion (tTw) must include all channels by which x can in- 
teract, including coannihilation with other particles, in 
which the number densities of both species are impor- 
tant. At some point in time the SUSY particle will no 
longer be able to remain in thermal equilibrium with its 
surroundings ( "freeze-out" ) and its co-moving number 
density will be nearly constant. We limit ourselves to 
models which obey a discrete symmetry, R-parity, which 
prevents decays (but not two-body scattering) of SUSY 
particles into standard model particles [1]. Any SUSY 
particle other than the LSP in existence at freeze-out 
will decay to the LSP state. Following the derivation in 
psl ] , the expression for the present relic density of a dark 
matter particle in terms of the cross section and freeze- 
out temperature is 

10~^^cm^ 

^dnih'^ ~ 0.3xfy/g^ — - — ^ — . (2) 

{av) 

Here Xf = m^/Tf where Tf is the freeze-out tempera- 
ture, and 5, is the number of effective relativistic degrees 
of freedom of all species contributing to annihilation. The 
important features here are that total annihilation cross 
section controls density, and that mass does not enter the 
equation except through a weak dependence in g^, and x / . 

B. Relic— Density Constraints 

The WMAP survey, when combined with recent obser- 
vations of large-scale structure, currently provides the 



best constraints on the quantity ridmh^, where fidm is 
the ratio of dark matter density to the critical den- 
sity pc = 1.88/i^ X 10~^^g cm^'^. For our analysis, we 
use the most recent third year WMAP data combined 
with Sloan Digital Sky Survey large-scale structure data 
[igl ] to arrive at the tightest constraint on relic density, 
ridmh'^ = 0.11lt°;°°^^, here h = 0.7091:°;°^^ being the 
Hubble parameter in units of 100 km/s/Mpc. We apply 
bounds equal to twice these la limits for the following 
analysis, 

0.096 < n^rnh'^ < 0.122. (3) 

These limits are sufficient to incorporate other recent 
measurements of fldmh^ [lOjIlii, which do not differ from 
our figure by more than ~ Icr. It should be noted that our 
results are not particularly sensitive to the relic density 
constraints, and there are larger sources of error involved 
in the calculation than the relatively minor variations in 
experimental determinations of fldmh^- 



C. The Supersymmetric Neutralino 

The dark matter candidate we address in this paper is 
the lightest supersymmetric neutralino, denoted x, in the 
context of the minimal supersymmetric standard model 
(MSSM). To express the neutralino mass states as a lin- 
ear combination of Higgsino, bino, and neutral wino par- 
ticle states, we diagonalize the mass matrix. 



Ml 


—rriz cos /3 sin 
rriz sin f3 sin 9^ 





M2 

TTiz COS (3 cos 9w 
—rriz sin (3 cos 9^ 



—rriz cos /3 sin 
•niz cos (3 cos 9y, 




mz sin (3 sin 0^ 
—niz sin /3 cos 6 





Here (3 is the ratio of vacuum expectation values be- 
tween the two Higgs doublets, is the mass of the , 
9^ is the weak mixing angle, and Mi, M2, p. are the 
U(l) and SU(2) gaugino and Higgsino mass parameters, 
respectively. The physical states then become 

X? = A,B + B,W^ + aH° + D,Hl (4) 

with Af + Bf + Cf + Df = 1. Here i = 1 to 4 is a parti- 
cle index that will be suppressed in cases where the LSP 
is being discussed. For our first models, we will be dis- 
cussing instances in which the LSP is entirely Higgsino- 
or wino-like, corresponding to the conditions C^+D^ « 1 
or « 1, respectively. The masses we are investigating 
are going to be strictly > 1 TeV, and the off-diagonal 
blocks in the neutralino mass matrix are < so mix- 
ings are not large, and therefore the mass of the LSP is 



tightly controlled by the least of the three mass param- 
eters. The Higgsino parameter p, and SU(2) parameter 
M2 also appears in the chargino mass matrix: 

f Nh V2m.sin/3\ 
^ \ V 2to«, cos/3 ^ J 

Again, the off diagonal parameters here are small com- 
pared to the mass scale of interest. Thus our Higgsino- 
and wino-type dark matter models come with a nearly- 
degenerate chargino (fermionic partner of charged Higgs 
and W bosons) built into the model at high energies. This 
chargino will account for a large degree of the total an- 
nihilation cross section for these two dark matter types. 
Incidentally, because we are interested in high-mass dark 
matter, > 1 TeV, we will not be addressing the possibil- 
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ity of bino-type dark matter in this paper. Because there 
is no degenerate chargino state in this case, the total an- 
nihilation cross sections and mass limits on bino dark 
matter tend to be much lower than the other two va- 
rieties, even with strong coannihilation from sfermionic 
particles. 



II. MASS LIMIT RESULTS 

A. Pure States 

To explore cases in which the LSP is a pure wino or 
Higgsino, all SUSY masses are set to high values (> 10 
TeV) except for either the Higgsino or SU(2) (wino) mass 
parameter. The other relevant particle which appears in 
this situation is a slightly heavier chargino, and in the 
Higgsino case a second nearly degenerate Higgsino-type 
neutralino. As the other particles are at a significantly 
higher mass scale, they will be thermally suppressed prior 
to dark matter freeze-out and will not affect relic density. 
Our results for this particular region show a monotoni- 
cally increasing relic density with increasing mass, with 
no dependence on the tan (i parameter. Our bounds from 
combined WMAP and SDSS for these cases are: 

0.99 < TO;^(TeV) < 1.12 (Higgsino) (5) 
2.10 < mx(TeV) < 2.38 (wino) (6) 

Our mass limit for a pure wino state is consistent with 
that mentioned by Profumo [s^], who quoted a function 
narr^h"^ = c(m^(TeV))'' with 0.0225 < c < 0.0255 and 
1.90 < 7 < 1.92. 



B. Coannihilation with a Sfermion 

In ordered to systematically test the effects of coanni- 
hilation with a Higgsino, we tested each sfermion mass 
parameter, originally set to high values, by shifting them 
down to the coannihilation region, to a mass rrica which 
is slightly larger than jjL. Beginning with the limits set by 
Higgsino-chargino coannihilation, we attempt to find re- 
gions where these processes allow a larger Higgsino mass 
by increasing the effective cross section for annihilation. 
Depending on the specific interaction strengths for pro- 
cesses involving this new particle, the relic density may 
be increased or decreased as rrica is brought lower, that 
is, coannihilation may have a positive or negative effect 
on the mass limit. In our low-energy effective supersym- 
metric model the sfermion masses are all free parame- 
ters. For our notation we write rriqi, mm, "ij^, m^^, and 
rriei, with i = 1,2,3 being the generation index, for the 
left-handed squark doublet, right-handed up and down 
singlets, left-handed slepton doublet, and right-handed 
slepton singlet, respectively. This gives 15 free mass pa- 
rameters to examine in this coannihilation calculation. 



Using MicrOMEGAs and confirming our results with 
DarkSUSY, we determine that only coannihilation with 
the third generation of squarks allows the Higgsino mass 
to be increased beyond the amounts in the previous sec- 
tion while conforming to experimental limits. This is 
shown graphically here in Figure [H which show regions 
in mca^fJ- space that fall within the 2 * Icr bounds on 
relic density. For these figures, all masses other than fi 
and rrica are set to 50 TeV. At this mass these particles 
are effectively removed from the early-universe Boltz- 
mann equation, as their number densities are exponen- 
tially suppressed prior to freeze-out. The ratio of Higgs 
vacuum expectation values tan P is set to 30, and the su- 
persymmetry breaking parameters Ai are set to 0. No 
other sfermion increases the maximum Higgsino mass 
when allowed to coannihilate. Also, the effect of coanni- 
hilation with two or more sfermions is not found to be 
cumulative in general. When m^s « fi, bringing down 
any other mass either increases ridmh^ or has a negliga- 
ble effect. We did not find any cases in which compound 
coannihilation with several sfermions had a substantially 
effect on the mass bound. 

The ratio of Higgs doublet vacuum expectation values 
tan (3 has a significant effect in the sfermion coannihi- 
lation region. Our previous analysis was done with a 
typical value of tan/3 = 30, but, with a higher value, 
the effects of coannihilation are increased. Third gener- 
ation sfermion right-hand singlet parameters u3 and d3 
can also have a slight enhancement effect in this regime. 
The optimal combination is found to be /x « niq^ = 
0.9mi3 = 0.9mj3. In this case, m-^ may grow as large 
as 1.8 TeV while remaining within experimental density 
bounds. This is the greatest mass which was possible un- 
der any combination of coannihilating sfermions in our 
Higgsino scenario. 

For a wino-type LSP, no sfermion coannihilation ar- 
rangement was found to have the effect of raising the 
mass bound; all sfermion coannihilation schemes cause 
the mass limit to decrease. Therefore the highest mass 
limit for this type of dark matter is that found in the 
previous section. 

C. Annihilation through a Massive Higgs 
Resonance 

Another mechanism through which the annihilation 
cross section of a neutralino LSP might be greatly en- 
hanced is via a heavy Higgs resonance, tua ~ 2rni5p. 
At the multi-TeV scale, the masses of the heavy CP- 
even Higgs, CP-odd Higgs, and charged Higgs are all 
approximately degenerate. 

niA ~ TTiH ~ Tn-H^ 

Under this arrangement, not only is the cross section of 
LSP annihilation enhanced by the resonance, but so are 
coannihilations between any nearly degenerate charginos 
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FIG. 1: Results for Higgsino-squark coannihilation. We show allowed regions in mj-/i space, where is the mass for 4 

different sfermion species shown here: q3, the left-handed doublet of 3rd-generation squarks in the lower left, u3 and d3, the 
right-handed singlet partner of the top and bottom quarks in the lower and upper right, and /3, the left-handed doublet of 
3rd-generation leptons in the upper left. The two lines indicate the upper and lower 2a bounds on relic density. The physical 
masses of the sfermions in each case are very nearly equal to our mass parameter, with corrections < 5 GeV that can be ignored. 
Thus, the contour lines must end at mj = /i, for the Ifiggsino to remain the LSP. The highest vertical point in the regions for 

u3 and q3 indicate the highest possible Higgsino mass. In the case of 13, the effect of coannihilation on the upper mass bound 
is strictly negative. The d3 parameter has more complex behavior, but the overall effect is negative. 



or next to lightest neutralinos. It is expected that the 
cross section will be dominated by the CP-odd Higgs 
channel, as the contribution from CP-even Higgses van- 
ish in the low velocity limit due to the requirement of CP 
conservation in the intermediate state Profumo [s^l 
analyzed this region of parameter space in minimal super- 
gravity (mSUGRA) and anomaly-mediated (mAMSB) 
SUSY breaking models with non-universal Higgs masses. 
He found upper LSP mass limit of approximately 5 and 
12 TeV for mSUGRA and mAMSB, respectively, utiUz- 



ing 2(7 WMAP bounds and tan/3 — 40. We have followed 
a similar program, taking the low-energy neutralino and 
chargino mass matrix inputs as free parameters, with the 
aforementioned constraint on the LSP mass, therefore ex- 
ploring over the vector space of neutralino and chargino 
mixings. For two multi-TeV neutralinos interacting at 
zero velocity, the cross section for annihilation through a 
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CP-odd Higgs is [nil, 

(--)-^Y.'^f\9Aff?^ '"'^f I, (7) 

where QAxx ^^"^ 9Aff are the vertex factors for the cou- 
phng between the Higgs and neutrahno and final-state 
fermion species /, respectively, and Ta is the Higgs 
width. The vertex factors appearing in a neutralino- 
Higgs or chargino-Higgs junction involve products of 
gaugino and Higgsino mixing factors, and are therefore 
sensitive to the exact choice of mass parameters. For two 
neutralino LSPs annihilating through a CP-odd Higgs, 
we have (53.] 

9Axx ^ {9B - 9' A){Csxn(i-D cos 13) (8) 
9AUU = '-^^ (9) 
9Aaa = (10) 

where the LSP composition is denoted by parameters 
A through D as in equation[4l Here 'u' refers to up-type 
quarks and neutrinos, and 'd' refers to down type quarks 
and charged leptons, and g and g' are the SU(2) and U(l) 
coupling constants. Because the gAxx vertex factor is de- 
termined by products of gaugino and Higgsino fractions, 
the largest factors and highest annihilation enhancements 
tend to occur when the neutralino is an even mixture of 
gau gino and Higgsino particle eigenstates. As derived in 
[34| . the maximum cross section for the CP-odd Higgs 
channel is at tan (3 ~ ^Jm^Jm^ ~ 6.4; our investigation 
of different values for tan (3 found the highest mass limits 
at approximately this point. 

Our results for this type of model can be seen in Figure 
[H We systematically calculated mass limits correspond- 
ing to the upper 2 * Icr experimental bound over a large 
number of neutralino mixtures. This was done by setting 
niA = 2/x « 2m;j, and scanning over values oi fi — Mi 
and /X — M2. These differences alter the particle state 
content of the LSP via the neutralino mass matrix. The 
final results are shown in terms of the fractional LSP 
composition. As tan/3 is also an independent, relevant 
parameter, we included this as a variable and searched 
from tan/3 = 2 to 50. The results from this scan are 
plotted as mass contours as a function of fractional neu- 
tralino composition. We show plots for 3 values of tan /3 
in Figured! 

III. DETECTION RATES 

For each of the model types thus discussed, we have 
calculated the resulting local gamma-ray flux from anni- 
hilations in the galactic center. For our halo model, we 
utilized the predictions of [l^l, together with their fidu- 
cial value for the density normalization at the maximum 
radius of the central black hole's sphere of influence, ~ 



2pc. This factor is highly uncertain and it is possible 
that it has been increased by baryonic infall. Dark mat- 
ter annihilation can also have an effect on the density 
profile of the innermost part of the cusp, changing the 
signal intensity by a multiplier refered to as the boost 
parameter j5J], an effect not considered here. A reader 
wishing to use a different normalization can simply raise 
or lower our flux measurements by the square of the nor- 
malization multiplier, or by the boost factor in the case 
of modifications to the inner halo profile. 

The total cross sections for all annihilations to contin- 
uous and line features in the spectrum were calculated 
using the DarkSUSY program. In order to show a typ- 
ical spectrum as seen by a detector, we have folded the 
distribution with a gaussian with an energy width of 15 
percent, as was done in [55]. This also allows a reason- 
able visual comparison to be made of the prominence of 
the line emission feature against the continuous output. 

In Figure [3] we show the gamma ray flux from 5 differ- 
ent models, corresponding to the scenarios we have dis- 
cussed. Clearly, with the fiducial halo normalization none 
of these annihilation models can account for a significant 
fraction of the flux observed by the H.E.S.S. telescopes 
and confirmed by MAGIC. However, the highly uncertain 
contribution to the central density from baryonic com- 
pression has not been taken into account in these results. 
As flux increases quadratically with particle density, a 
rather modest compression factor of 10 would increase 
flux by 2 orders of magnitude, enough to bring our more 
strongly annihilating models to the levels observed by 
these atmospheric Cherenkov telescopes. However, the 
H.E.S.S. data also maintains an approximate power-law 
proflle for 2 logarithmic decades, something that none of 
our models can reproduce even with a carefully adjusted 
density normalization. Even the most massive particles 
we found to be capable of satisfying relic density con- 
straints exhibit a roll-off behavior that is not in the ob- 
served spectrum. 

Two plots for the annihilation through a heavy Higgs 
resonance are shown. We have chosen a mass of 20 TeV 
here, the heaviest neutralino for which a cross section 
could be computed without resorting to extrapolation in 
certain DarkSUSY routines. Because mean halo veloci- 
ties are much lower than those at freeze-out, annihilation 
cross sections are more sensitive in the halo to the exact 
relation between ruA and as there is little smearing 
out of the center-of-mass energy due to thermal veloc- 
ities. Thus models that yield similar relic densities can 
have very different halo annihilation cross sections. To il- 
lustrate this we have displayed both an optimized model 
(orange) with mA = 2m-^ and a more typical model (pur- 
ple) where the relation is only approximate. 



IV. DISCUSSION 

Higgsino- and wino-type dark matter annihilate much 
more strongly than bino dark matter at a given mass. 
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FIG. 2: Resonant heavy Higgs annihilation mass hmit plots for tan/3= 2 (top), 10 (middle), and 50 (bottom). The axes are 
the fractional wino and bino states of the LSP, with the Higgsino fraction being the remainder. The contours show the upper 
mass bounds in TeV. 
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FIG. 3: In the upper plot, we summarize our findings by showing the resulting local gamma-ray flux from the galactic center in 
several annihilation scenarios using the halo model of with fiducial normalization (no baryonic compression), and compare 
to the latest observations of the H.E.S.S. experiment (black data points, [30])- The dashed lines show the true continuous 
distribution, while the solid lines show the total (continuous plus discrete) emission spectra as seen by a detector with an 
energy resolution of 15 percent. The blue line is a 1 TeV Higgsino, coannihilating with a nearly degenerate chargino and second 
Higgsino. The red line shows the same model with coannihilation from a 3rd generation squark, at a mass of 1.8 TeV. The 
green line is a 2.4 TeV wino. The purple and orange lines are both a mixed type neutralino annihilating through a heavy 
Higgs resonance. The orange model has been optimized by fine tuning of the resonance, so that the cross section and resulting 
fiux are maximized, while the purple line shows a more typical model. The lower plot demonstrates an attempt to fit a Higgs 
resonance model to the H.E.S.S. data. A factor 10 density boost is applied, resulting in a 10^ increase in flux above the fiducial 
value. 



and are capable of satisfying relic density constraint with 
masses at the TeV scale, something that is not possible 
with a bino particle, even including sfermion coannihila- 
tions. Primarily Higgsino-type dark matter can arise, for 
example, in the 'focus point' region of minimal supergrav- 
ity (mSUGRA) [5^. In this model, the supersymmetric 
scalar quarks and leptons are set to high masses by a 
single GUT scale parameter, and the neutralino becomes 
a Higgsino-bino mixture. As this scalar mass term is 
increased, the neutralino LSP becomes more Higgsino- 
like, and features an increasing cross section [Ij, ^t\ . 



Since the highest mass LSPs in the focus point region are 
nearly pure Higgsino, we have constrained this and other 
Higgsino-type models by considering this limiting case. 
Dark matter which is predominantly wino-type appears 
in minimal anomaly-mediated supersymmetry breaking 
[58| . MicrOMEGAs provides information on the annihi- 
lation channels relevant to the total LSP cross section in 
a particular model. For the case of Higgsino-chargino 
coannihilation, a large number of channels provide mod- 
est contributions to the total cross section. Chargino and 
neutralino annihilation to quarks are the most impor- 
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tant processes, with annihilations to leptons and gauge 
bosons, and double chargino annihilation to quarks being 
the other relevant channels. In the >TeV mass range of 
interest, even the heaviest standard model particles are 
essentially massless, and there is little difference in avail- 
able phase space between different generations, hence lit- 
tle variation in cross section. Not surprisingly, there is 
also relatively little variation in channel contributions as 
the mass scale and resulting relic density are altered. 

The introduction of a coannihilating squark opens the 
possibility of tree-level annihilations to gluons. Annihila- 
tions to gluons, as well as coannihilations between Higgsi- 
nos and squarks to quarks and gauge bosons, are the ma- 
jor new channels available. The decreased mass of a given 
sfermion increases the t-channel amplitude for scattering 
from neutralinos to fermions, but only for that particular 
flavor of the sfermion. These same vertices explain why 
the third generation of squarks are unique as the only ef- 
fective set of coannihilating partners. The appearance of 
the corresponding quark mass terms limits these channels 
to cases involving heavy quark masses. The quantity (3 
appears in several vertex factors involving Higgs and Hig- 
gsino iterations with matter [sl] . Processes with a factor 
of tan/3 or (sin/3)~^ in the amplitude include chargino 
and neutralino annihilation with a sfermion to a gauge 
boson and fermion, and these are primarily responsible 
for the decrease in relic density in scenarios with squark 
coannihilation with increasing tan j3. 

An effect which is relevant for Higgsino and wino mod- 
els is the Sommerfeld enhancement, appearing in the 
context of weakly-interacting non-relativistic particles 
with a mass much greater than the W-boson. This non- 
perturbative effect has been shown to decrease relic abun- 
dance by as much as 50% for a wino-type LSP and 10% 
for a Higgsino-type [s^ . These authors study wino dark 
matter as an example and find upper WMAP bounds of 
2.7 to 3 TeV. Neither DarkSUSY nor MicrOMEGAs ac- 
counts for this effect; in addition to these considerably 
higher new bounds for wino dark matter, our bounds in- 
volving pure Higgsino dark matter could be enhanced by 
a slight amount. 

When annihilation through a heavy Higgs s-channel 
resonance is considered, allowed masses can go into the 
tens of TeV. The examination of this scenario was done 
using the DarkSUSY software. As expected, the neu- 
tralinos which had the largest cross section were approx- 
imately even mixtures of Higgsino and gaugino states. 
From Figure [2l it is clear that there is little change in 
the general topography of the relation between neutralino 
mixture and mass bounds with changing tan /3. The max- 
imum mass limits do change as a function of tan (3, rising 
a small amount from 32 TeV at tan/3 = 2 to about 34 
TeV at tan/3 = 5 to 8, and then decreasing from that 
point down to 18 TeV at tan /3 = 50. 

Our results with MicrOMEGAS were tested against 
the DarkSUSY code, and the programs were found to 
generally be in agreement over the parameter space of 
interest, except in the case of annihilation through a 



Higgs resonance. For our calculations involving pure Hig- 
gsino and wino states, the difference in relic density cal- 
culations were no larger than 2.5 percent, and in cer- 
tain cases where sfermion coannihilation was considered 
it was no greater than 9 percent. It was noted during 
while investigating sfermion coannihilation that the two 
codes produced highly disparate results in certain situa- 
tions. These problems were determined to be a error in 
the DarkSUSY software that only appears at mass scales 
higher than we have considered here, and did not appear 
to be an issue for our results in the coannihilation region. 
For the Higgs resonance models, there was a significant 
difference in the predictions of the two codes, sometimes 
by as large as a factor of 2. The DarkSUSY results, which 
we have presented, tended to output higher mass limits 
than MicrOMEGAs. While this does mean that our re- 
sults in this area should be taken only as approximate, 
our conclusion that this scenario is unlikely and cannot 
explain current gamma ray observations is not altered. 

While there is no concrete upper bound on the scale 
of supersymmetry breaking, a mass well into the TeV 
range is certainly disfavored by constraints from gauge 
coupling unification [63,[6l|. Another 'absolute' bound 
comes from partial wave unitarity, which provides an up- 
per limit on the mass that any thermally produced dark 
matter particle can have, by placing a constraint on the 
cross section in equation [T]. This bound is applicable 
as long as the annihilation cross section arises primarily 
from s-wave terms. The mass limit set by unitarity is 
^dmh"^ > 1.7* 10-^ y^{mdm/TeVf [63|, which leads to 
a maximum relic density of about 120 TeV. While this 
mass is well above that of any of the MSSM models we 
examined, it may become important when considering 
thermally-produced heavy dark matter candidates from 
other particle physics extensions to the standard model. 
However, it should be noted that the unitarity bound can 
be violated in the case of a strong resonance, in which 
the assumption of s-wave dominance breaks down [6^. 
Another case in which the unitarity bound is not applica- 
ble comes from possible non-perturbative factors in the 
annihilation cross section which could affect heavy (>500 
GeV) Higgsino-type neutralinos [U- These factors ap- 
pear only at low velocity and would not affect the physics 
of the dark matter during freeze-out but could thus af- 
fect halo interactions, greatly increasing flux levels from 
annihilations. 

The mass limits we have set in this paper apply only 
to neutralino dark matter in the MSSM model which at- 
tains a relic density through thermal freeze-out in a stan- 
dard cosmology. We can make no claims about cases in 
which the dark matter is produced through non-thermal 
processes. These could include scenarios in which the 
dark matter is produced non-thermally, possibly by a 
late-decaying scalar field or one in which entropy 
is produced after freeze-out. This latter case could hap- 
pen for a variety of reasons ([^, for review) and would 
have the effect of violating the standard assumption of 
constant comoving entropy density, which would reduce 
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relic density. The unitarity bound would not apply in this 
situation, as a very massive dark matter particle with or- 
dinary cross section could still attain the correct density 
today. 

V. SUMMARY 

We have determined the masses of pure Higgsino and 
wino-type thermally-produced dark matter which are 
consistent with the latest density constraints on dark 
matter, defined here as twice the la boimd determined 
by combined SDSS and WMAP-3 data. In the absence 
of any coannihilation processes with scalar fermions, the 
suitable mass range is found to be between 0.99 and 1.12 
TeV for a pure Higgsino and 2.10 and 2.38 for a pure 
wino state. Coannihilation with partners of the 3rd- 
generation quarks is found to increase this limit modestly 
for Higgsino type dark matter to an upper limit of about 
1.80 TeV, with fine tuning in the mass parameters and 
tan/3, but no coannihilation model can increase the mass 
limit for a wino type particle. Allowing the dark matter 
to exist as a bino pure or mixed state tends to sharply 
decrease mass limits, and bino mass limits were always 
found to be in the sub-TeV range. The other class of 



models which wc examined utilized annihilation of the 
LSP through a heavy Higgs resonance. Viable models 
with LSP masses as high as 34 TeV were found, though 
these scenarios are sensitive to the both the neutralino 
mixture and the resonance condition niA = 2m j^, and 
are therefore dependent on fine-tuning. A computation 
of the VHE gamma ray spectrum which could be ob- 
served with an atmospheric Cherenkov telescope showed 
that even the largest masses wc found arc not adequate 
for fitting the observed H.E.S.S spectrum. This observed 
event rate is also considerably higher than our predic- 
tions, although this is not necessarily a problem because 
of uncertain normalization of the dark matter profile in 
its innermost regions. 
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